Formation of guided spin-wave bullets in ferrimagnetic film stripes 
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The formation of quasi-2D nonlinear spin- wave eigenmodes in longitudinally magnetized stripes of 
a ferrimagnetic film, so-called guided spin-wave bullets, was experimentally observed by using time- 
and space-resolved Brillouin light scattering spectroscopy and confirmed by numerical simulation. 
They represent stable spin-wave packets propagating along a waveguide structure, for which both 
transversal instability and interaction with the side edges of the waveguide are important. The 
experiments and the numerical simulation of the evolution of the spin- wave excitations show that the 
shape of the formed packets and their behavior are strongly influenced by the confinement conditions. 
The discovery of these modes demonstrates the existence of quasi-stable nonlinear solutions in the 
transition regime between one-dimensional and two-dimensional wave packet propagation. 
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Stable two-dimensional localized nonlinear spin-wave 
excitations, so-called spin-wave bullets, have been pre- 
viously observed in thin ferrimagnetic films of yttrium- 
iron-garnet (YIG) magnetized along the propagation di- 
rection 0, 2, 3]. These films were practically unbounded 
in both in-plane directions compared to the transver- 
sal size of the spin-wave packets and the wavelength of 
the carrier spin wave. In contrary, in a one-dimensional 
waveguide structure, where the width is comparable to or 
smaller than the spin-wave wavelength, only quasi one- 
dimensional nonlinear spin-wave objects were observed, 
which are spin- wave envelope solitons. 0, Sj- Both for 
solitons and bullets linear pulse spreading in the direc- 
tion of propagation (so called longitudinal direction) due 
to wave dispersion is compensated by longitudinal non- 
linear compression. As for the transverse in-plane direc- 
tion, solitons are meant to have a stable transverse dis- 
tribution of their dynamic magnetization coinciding with 
the profile of the lowest linear spin-wave eigenmode of 
the waveguide. On the contrary, bullets show transverse 
nonlinear instability of attractive type which overcom- 
pensates transverse diffraction broadening of the wave 
packet. The nonlinear compression would lead to a wave 
packet collapse if the medium is lossless. Weak mag- 
netic losses in a real magnetic film ensure a fine balance 
of nonlinear narrowing of the packet and of its diffrac- 
tion spreading for some distance of propagation [l[ . This 
results in a quasi-2D spatially localized bell-shaped wave- 
form which is stable during the lifetime of the bullet. 

Here we report on the experimental observation of a 
stable spin-wave packet propagating along a waveguide 
structure, for which both transversal instability and in- 
teraction with the side edges of the film waveguide are 
crucial. The structure can be considered as a transitional 
case between the 2D case of a continuous film and the 
quasi-lD case of a narrow stripe. We show that in this 
case the nonlinear wave dynamics is distinctively differ- 
ent from both the soliton and the bullet cases, and can be 



considered as efficient coherent nonlinear mixing of spin- 
wave eigenmodes of the waveguide. Our theory shows 
that this mixing is due to a specific interaction - the 
pseudo-linear generation of higher order modes by the 
fundamental one. 

The experiment was carried out using a longitudinally 
magnetized long YIG film stripe of 2.5 mm width and 
7^tm thickness. The magnetizing field was 1831 Oe. The 
spin waves were excited by a rf magnetic field created 
with a microstrip antenna of 25 /xm width placed across 
the stripe and driven by input microwave current pulses 
20 ns in duration at a carrier frequency of 7.125 GHz. 
The spatio-temporal behavior of the traveling spin-wave 
packets was investigated by means of space- and time- 
resolved Brillouin light scattering spectroscopy Q. 

The results of our measurements are demonstrated in 
Fig. [T] The panels show the spatial distribution of the 
intensity of spin-wave packets at different times of their 
propagation from the left to the right along the stripe 
waveguide. The left vertical set of diagrams corresponds 
to the linear case. The input microwave power for this 
set is 20 mW. The right set is for nonlinear propaga- 
tion. These data were collected applying an input driving 
power of 376 mW. 

Differences between these two cases are clearly seen. 
The linear spin-wave packet is characterized by a trans- 
verse profile very similar to one half of the period of the 
sine function, while the cross section of the nonlinear 
packet has a pronounced bell-like shape. Furthermore, 
the intensity of the linear packet decays monotonically 
with time. That is obviously due to magnetic damping 
in the film. The intensity of the nonlinear packet ini- 
tially increases because of its strong transverse compres- 
sion (see the second diagram from the top in Fig.[T]), then 
gradually decays as the damping overcomes the nonlinear 
compression for larger propagation times. 

Formation of a stable nonlinear spin-wave object is also 
evidenced in Fig. [2] which shows the measured peak in- 
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FIG. 1: Observation of linear propagation and of formation 
of the guided spin-wave bullet in an YIG film stripe waveg- 
uide using time- and space-resolved Brillouin light scattering 
spectroscopy. For parameters see main text. 



tensity of the nonlinear wave packet and its width (i.e. 
packet size in the transverse direction) as a function of 
propagation path. At the beginning of the propagation 
path one clearly sees an increase of the peak intensity 
on top of the exponential decay. The increase in peak 
intensity is due to transverse compression of the packet, 
evidenced in the same figure by a decrease of the packet 
width. Down the propagation path the width stabilizes, 
and in the range of distances z=l-4mm it does practi- 
cally not change, apart from small superimposed oscilla- 
tions which will be discussed below. 

Both figures provide a clear proof of development of 
transverse instability and bullet formation. Interestingly, 
the 2D bell-like shape survives even at the end of the 
propagation path (z >4mm) when the packet intensity 
has decreased more than ten times and the nonlinearity 
contribution to the spin- wave dynamics has been consid- 
erably diminished. 

Importantly, in Fig. [5] one observes superimposed oscil- 
latory variations in the intensity and in the packet width, 
both in the quasi-linear regime of the packet propagation 
(z > 4 mm) and in the highly nonlinear regime z < 4 mm. 
In the entire range the oscillations in intensity and in 
width are in anti-phase to each other. This effect has not 
been observed in case of conventional bullet formation [l[ . 
On the contrary, a similar picture is usually observed for 
linear guided spin waves in narrow film waveguides 0, S] j 
where it is explained as a beat of phase-correlated linear 
waveguide width modes propagating at the same carrier 
frequency. This effect evidences the importance of the 
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FIG. 2: Measured peak intensity and width for the nonlinear 
wave packet shown in the right panels of Fig. [1] as a function 
of propagation path. The dashed vertical lines show positions 
of local minima of intensity and of the local maxima of the 
packet width. 



influence of confinement on the nonlinear evolution of 
the packet transverse profile and suggests that interac- 
tion of the linear eigenwaves of the waveguide - the so- 
called width modes - underlies nonlinear wave dynamics. 
Therefore we term this nonlinear object a "guided spin- 
wave bullet" . 

The theoretical description of the observed phenomena 
is based on concepts developed in Refs. 0| and [Io[. The 
ingredients of the model are (i) effective dipole pinning 
of magnetization, which results in a tri-linear interaction 
as the initial interaction for development of the trans- 
verse instability of the wave packet profile, (ii) the width- 
mode group velocity matching, and (iii) the nonlinear 
extension of spectrum of width modes. In Ref. 0], we 
theoretically studied linear propagating eigenwaves of a 
magnetic stripe. The eigenwaves represent guided modes 
with discrete transverse wavenumbers. For a stripe of 
rectangular cross-section the modes are characterized by 
a standing-wave type dynamic magnetization distribu- 
tion across the stripe cross-section and by a monochro- 
matic propagating wave with the longitudinal wavenum- 
ber k z along the stripe. Importantly, for stripes with a 
large ratio p of width to thickness the thickness distri- 
bution of the dynamic magnetization is practically ho- 
mogeneous, whereas in the direction of the stripe width 
the standing spin waves possess considerably decreased 
amplitudes at the edges due to dynamic demagnetiza- 
11 1 . Previous calculations (see Fig. 3 in [§]) 



tion effects 

show that the assumption of a totally pinned magneti- 
zation at the stripe edges results in good approximation 
for the transverse profiles of the propagating eigenmodes. 
Adopting the assumption of totally pinned edge spins, 
the transverse profile of the dynamic magnetization is de- 
scribed by an integer number of half-periods n of the sine 
function with the lowest transverse wavenumber k y = k n 
being fc n =i = 1 • tt/w, where w is the stripe width. For 
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the aspect ratio of our waveguide, p = 357, this works 
with good accuracy. 

As one sees from Fig. [TJ this conclusion is in a good 
agreement with the packet profiles measured in the lin- 
ear regime. The nonlinear dynamics is now described 
by a theory which is analogous to the one developed in 
Ref. [13]. We assume that the dynamic pinning of mag- 
netization is conserved in the weakly nonlinear regime. 
Then the transverse evolution of the nonlinear packet can 
be considered as interaction of linear eigenmodes which 
are pinned at the stripe edges. Indeed, in Fig. [T] one 
clearly sees that in the nonlinear regime the dynamic 
magnetization at the stripe edges practically vanishes. 
Then the description for the propagating modes results 
in an evolutional equation for the spin-wave precession 
angle <f> which reads: 

id/dtF n ,k[<i>(y,z)] + (u n ,k +irj- u)F n j<.[<l){y,z)] (1) 
+TF n>k [\ <j>{y,z) | 2 tj>(y,z)] = f n ,k,t- 

In this expression k = k z , uj n ,k is the eigenfrcquency of 
the n-th guided mode for the longitudinal wavenumber k, 
r\ is the relaxation frequency for the film, and u> is the car- 
rier frequency of the microwave signal f(y, z, t) exp(iwi) 
applied at the vicinity of z—0 which excites the input 
spin- wave packet, and T is the nonlinear coefficient anal- 
ogous to the nonlinear coefficient of the spin-wave version 
of the Nonlinear Schrodinger Equation [12j . The opera- 
tion F n> k denotes a 2D Fourier transform. It is the dis- 
crete sine transform with the basis functions sin(fc n y) in 
the direction of stripe width y, and is the Fourier integral 
over continuous wavenumbers k with the basis functions 
exp(ikz) along the stripe. 

The expression Eq. (1) can be easily transformed into 
a system of dynamic equations for amplitudes of guided 
modes <p n (z,t) — F n [<p(y, z,t)] coupled by the four-wave 
nonlinear interaction. The analysis of the system shows 
that the formation of the two-dimensional waveform can 
be considered as an extension of the spectrum of the 
width modes. The partial waveforms carried by the in- 
dividual width modes have the same carrier frequencies 
equal to that of the external excitation signal to and the 
carrier wave numbers which satisfy the dispersion rela- 
tions for the modes io n ^ = lu. In the linear regime all the 
modes are independent. In the nonlinear (high ampli- 
tude) regime the width modes become mutually coupled 
which ensures intermodal coherent energy transfer. 

The efficiency of mode coupling in the pulse regime 
depends on two major factors: the mode group velocity 
matching and the type of nonlinear interaction. The ge- 
ometry of a relatively wide stripe is very favorable for 
having maximum contributions from both. 

Let us first discuss the type of nonlinear mode inter- 
action. The spin-wave packet immediately after having 
been launched into the stripe is carried by the lowest 
(fundamental) width mode (n = 1) <j>\(z,t). Therefore 



it is necessary to consider the nonlinear interaction of 
higher-order width modes with this particular mode. 

The interactions of the fundamental width mode with 
all even modes is not important for symmetry reasons. 
The nonlinear interaction of modes of the same type of 
symmetry is described by the parametric term as well 
as by an additional pseudo-linear (tri-linear) term. The 
parametric interaction is of convective instability type 
and is the same parametric instability which triggers for- 
mation of conventional bullets in continuous films. This 
conventional process is described by a pair of complex- 
conjugated equations with a parametric term propor- 
tional to the square of amplitude of the pumping wave 
(see e.g. Eqs. (4) and (5) in [13]). In this interaction 
the packet carried by the fundamental mode plays the 
role of the pumping wave. Its energy is transferred to 
the partial waveforms carried by the higher-order width 
modes. There is a threshold associated with this para- 
metric process due to natural damping in the medium. In 
the waveguide structure the threshold is of the same order 
of magnitude as the modulation instability threshold in 
continuous films (see e.g. Eq. (10) in Ref. Q). Further- 
more, an initial perturbation in the form of non- vanishing 
amplitude of a higher-order mode is needed to start the 
parametric amplification. This perturbation usually is 
provided by thermal excitation. Therefore an amplified 
higher-order mode, which is group velocity matched with 
the pumping wave, needs a large distance of propagation 
in order to reach a noticeable level. The energy of the 
pumping wave decreases down the propagation path due 
to losses in the medium. If the parametric amplification 
gain is small because of small supercriticality the higher- 
order mode cannot reach an amplitude comparable with 
that of the fundamental mode before the amplitude of the 
pumping wave packet carried by the fundamental mode 
falls below the threshold of parametric instability and the 
gain ceases. 

What distinguishes the confined waveguide geometry 
is that the nonlinear mixing starts as a pseudo-linear 
(tri-linear) interaction of the fundamental with the next- 
order symmetric mode which is the third mode <f>3(z,t): 
dfa/dt + v 3 d(j)3/dz + iuj^fa — Siz(z, t). Here v% is the 
group velocity of the third mode, wf 1 is its nonlinear 
frequency shift, and S13 is the tri-linear inhomogeneous 
term. This term has the form of a linear source of exci- 
tation with amplitude proportional to |</>i(z,i)| 3 moving 
with the group velocity of the fundamental mode. The 
presence of this pseudo-linear interaction at the early 
stage of the bullet formation is entirely due to the ef- 
fective dipolar pinning of the magnetization at the stripe 
edges. If the edge spins were unpinned, the interaction 
of all the width modes would be purely parametric. The 
pseudo-linear excitation introduced by this term into the 
dynamic equations is a threshold- free process, so it works 
even below the threshold of parametric instability. Fur- 
thermore, in contrast to the parametric process this pro- 



cess does not require non-zero initial amplitude of the 
amplified waveform to start the amplitude growth. This 
mechanism ensures rapid growth of the symmetric n = 3 
mode driven by the n — 1 mode up to the level where the 
parametric mechanism starts to work efficiently. After 
that the fundamental mode jointly with the n = 3 mode 
is capable to rapidly generate a large set of modes with 
yet higher odd numbers n through both pseudo-linear 
and parametric mechanisms. 

Let us now discuss the importance of the mode group 
velocity matching. Our theory shows that the efficiency 
of both nonlinear interaction mechanisms (parametric 
and tri-linear) strongly depends on the group velocity 
difference of interacting modes and the initial length of 
the nonlinear packet. In wider stripes the group veloc- 
ities of the width modes are closer to each other. As 
the nonlinearity is of attractive type, nonlinear correc- 
tions to the group velocities partially compensate for the 
group velocity mismatch. As a result the nonlinearly gen- 
erated higher-order partial waveforms remain for some 
time within the pump packet. Again, due to the attrac- 
tive character of the nonlinearity the modes have initial 
phases such as their transverse profiles are summed up 
constructively in the middle of the stripe width and a 
bullet-like total wave packet is formed as confirmed by 
direct numerical solution of Eq.(I) shown in Fig. [3j 

In the narrower stripes (I mm in width as in Fig. [31 
panel 2) the mode group velocity difference is larger and 
cannot be compensated by the nonlinear corrections to 
the group velocities for the same initial intensity of the 
wave packet. Thus the group velocity matching is not 
ensured. As a result the nonlinearly generated higher- 
order modes leave the area of the interaction before they 
reach significant amplitudes. For the same length of the 
initial packet <f>\ (t — 0, z) the extension of the spectrum 
of width modes does not occur. The waveform <fii (t, z) 
remains unaffected by the interactions with the higher- 
order modes. It undergoes only a nonlinear longitudi- 
nal compression and forms a quasi-ID wave packet - the 
spin-wave envelope soliton - which has a stable sine-like 
profile in the y-direction. 

The excellent agreement of the simulation results with 
the experimental data shown in Fig. [3] provides evidence 
for the validity of the developed theory. 

In conclusion, the formation of quasi-2D nonlinear lo- 
calized wave packets - guided spin-wave bullets - was 
studied in spin-wave waveguides. Our experimental and 
theoretical investigations show that formation of these 
stable nonlinear objects is strongly affected by the trans- 
verse confinement of the medium. A specific magneto- 
static effect - the effective dipole pinning of the mag- 
netization at the edges of the stripe, the width-mode 
group velocity matching of different discrete waveguide 
modes, and the extension of the width mode spectrum 
due to nonlinear mode-mode energy transfer are essen- 
tial for the nonlinear evolution of the initial spin-wave 
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FIG. 3: Lateral shapes of the nonlinear SW packets, f and 2 
- theoretical results calculated for the ferrite stripes of width 
of 2.5 mm and f mm, respectively. 3 and 4 - experimental 
profiles observed in YIG waveguides of width of 2.5 mm and 
f mm, respectively. 1 and 3: bullets. 2 and 4: solitons. 



excitation. Both the experimentally detected properties 
of the evolution and the theoretically revealed mechanism 
of formation show that the observed nonlinear wave pack- 
ets can be treated as "guided spin-wave bullets" which 
are specific for laterally confined magnetic films. 
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